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ABSTRACT 

We examine all possible stationary, optically thick, geometrically thin ac- 
cretion disc models relevant for active galactic nuclei (AGN) and identify the 
physical regimes in which they are stable against the thermal-viscous hydrogen 
ionization instability. Self-gravity and irradiation effects are included. We find 
that most if not all AGN discs are unstable. Observed AGN therefore represent 
the outburst state, although some or all quasars could constitute a steady pop- 
ulation having markedly higher fuelling rates than other AGN. It has important 
implications for the AGN mass supply and for the presence of supermassive black 
holes in nearby spirals. 

Subject headings: Subject headings: accretion, accretion discs — galaxies: active 
— instabilities 
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1. INTRODUCTION 

Accretion discs are a nearly ubiquitous feature of close binary systems, and their pres- 
ence is widely invoked in models of active galactic nuclei (AGN). A major feature of the 
discs in binaries is the thermal-viscous instability driven by hydrogen ionization (Meyer & 
Meyer-Hofmeister 1982, Smak 1982). It is now commonly accepted that this instability 
drives outbursts in cataclysmic variables and soft X-ray transients. Lin & Shields (1986) 
showed by a local stability analysis that this instability can also operate in accretion discs 
thought to be present around supermassive black holes in AGN. They concluded that these 
discs were unstable at radii 10 15 — 10 16 cm), where the surface temperature is several 
thousand degrees. The expected characteristic time scale for this instability is 10 4 — 10 7 
years. 

Because of its generic nature, the ionization instability plays a dominant role in char- 
acterizing the observed behaviour of the host systems. In the binary context, attempts to 
understand the precise conditions (mass of accreting object, accretion rate) under which it 
occurs have been at least partially successful (e.g. Smak, 1982, van Paradijs 1996; King, 
Kolb, & Burderi, 1996; King, Kolb & Szuszkiewicz, 1997, and references therein). These 
studies show that self-irradiation of the disc by the central X-ray source has a determining 
effect on the disc stability if the accreting object is compact, as in soft X-ray transients (see 
below). Delineating the stable and unstable disc regions is equally important for AGN. If 
the instability is present in AGN discs the suppression of central accretion in the quiescent 
state means that we can identify only the outburst states of unstable systems as AGN. Two 
important consequences follow: (1) quiescent AGN must appear as quite normal galaxies, 
and (2) the average mass fuelling rate in many if not all AGN is much lower than implied 
by their current luminosities. This in turn limits the masses that their central black holes 
are expected to reach. 
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The idea of intermittent activity in AGN was already suggested by Shields & Wheeler 
(1978). They noticed that the fuelling problem could be solved if active nuclei store mass 
during quiescence and this mass then feeds the hole for a shorter period of intense activity. 
The thermal-viscous hydrogen ionization instability found to operate in AGN accretion 
discs (Lin & Shields, 1986) is capable of triggering such behaviour. Clarke & Shields (1989), 
Mineshige & Shields (1990), Cannizzo & Reiff, 1992, Cannizzo (1992) studied the full range 
of black hole masses and accretion rates in order to determine the observational consequence 
of the instability for the AGN population. Siemiginowska & Elvis (1997) attempted to 
reproduce the observed luminosity function, assuming that this mechanism operates in all 
AGN. 

Our aim here is to decide if the ionization instability still operates in AGN when ir- 
radiation effects are included. As we have seen irradiation is central to the discussion of 
disc stability in soft X-ray transients. Further, irradiation is often thought to dominate the 
disc emission (e.g. Collin-Souffrin, 1994). For both reasons it is vital to include it in any 
attempt to decide the disc stability. The actual form of the instability when irradiation is 
included is outside the scope of our paper. Siemiginowska, Czerny & Kostyunin (1996) have 
performed studies for particular black hole masses and accretion rates, with assumed forms 
of irradiation. 

There is a simple criterion for the instability to appear: the disc must contain regions 
with effective temperature T eS close to the value T H at which hydrogen is ionized. In practice 
T H depends on the density and may be quite different in different environments; we shall 
consider a range of values in this paper. However the criterion is not easy to use in this form, 
as one does not in general know the radial distribution of the accretion rate, and thus the run 
of T e g, in a time-varying disc. Accordingly one usually uses the criterion in an indirect form: 
a disc with a given constant accretion rate M is self-consistently steady only if T eS > T H 
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throughout it. If the criterion fails we may expect outbursts, although the precise nature of 
these will depend for example on the detailed behaviour of the disc viscosity. 

This version of the stability criterion is easy to apply. Since T e g always decreases with 
disc radius R in a steady disc, the condition is most stringent at the outer disc radius R out , 
so we need apply it only there. If the disc's only source of energy is local viscous dissipation 
we have 

4 _3GMM 

[Teff(jR)] " (1) 

(e.g. Frank, King & Raine (1992); all symbols are explained after equations (2-9)), and 
the criterion is simply T e s(R out ) > Tr. In a binary system we can estimate R out with 
reasonable accuracy as 70% of the Roche lobe radius of the accreting star, and the problem 
is now well determined. Using this approach, Smak (1982) successfully divided outbursting 
cataclysmic variables (dwarf novae) from the persistent systems (novalikes). The extension 
to low-mass X-ray binaries is complicated by the fact that the dominant heat source for the 
disc is not local viscous dissipation (equation (II])), but irradiation by the central X-rays. 
The instability is similarly suppressed if the disc surface temperature given by irradiation 
exceeds Th (Tuchman, Mineshige & Wheeler, 1990). Provided that due account is taken 
of this, one can again successfully divide the outbursting systems (soft X-ray transients) 
from the persistent systems (van Paradijs, 1996; King, Kolb & Szuszkiewicz, 1997). The 
key feature, as in the unirradiated case, is that the edge temperature of the disc can be 
simply expressed in terms of M and R on t, without any need to solve for the full internal disc 
structure. In both the CV and LMXB cases there are important consequences for the study 
of the binary evolution (e.g. King, Kolb & Burderi, 1996), which gives a connection between 
M,M and R out . 

The extension of this approach to AGN is more complicated; here the outer edge of the 
disc is no longer determined by the simple Roche lobe condition which holds in binaries, but 



by the requirement that the disc becomes locally self-gravitating (see equation (|13D below). 
This condition requires a knowledge of the disc density at the outer edge, so we are now 
required to solve the full global structure of the steady disc to find -R ou t- Thus we examine 
all possible, stationary, optically thick, geometrically thin disc models relevant for AGN. If 
these correspond to stable states, AGN discs may be globally steady, and require fuelling at 
the currently inferred central accretion rates. If not, they will be the outburst states, and 
the required fuelling rates will be lower than the current central accretion rate. 



2. GLOBAL DISC STRUCTURE 

As explained above, to apply the stability criterion we need the self-gravity radius, 
and thus the global structure of steady discs. We assume these to be optically thick and 
geometrically thin. If we exclude the region close to the central object and consider total 
disc luminosities L ^ 0.2Le (£e being the Eddington luminosity) these approximations are 
justified, and we may parametrize all possible disc structures by M, M and the viscosity 
parameter a (see below). Our approach exploits only hot radiative regime of the disc, for 
which vertically average structure is a good approximation. The algebraic system describing 
a Shakura-Sunyaev disc is (see e.g. Frank, King, & Raine, 1992, p. 77): 

p = EH' 1 (2) 



H = CsR^G-^M- 1 ? 



(3) 



4aT 4 3GMM p „ 



t = Ere (7) 



E = — Mz/- 1 f 
3tt 



2 

z/ = -aHc s for a— discs (9) 



2 / P \ 
z/ = -aHc s \ f° r discs. 



Here G is the gravitational constant, R is the radial distance from the central object (assumed 
not to be extremely close to the inner radius of the disc), E is the surface density of the disc, 
ft — ftabs + ft es is the opacity (re a bs being the part due to (mainly) free-free and free-bound 
processes, and k cs being the part due to electron scattering), H is the disc half-thickness, c s 
is the sound speed, p is the density of the disc, T is its central temperature related to the 
effective temperature by T 4 = (3r/4)T c 4 ff , Q lvv is the flux absorbed from the external radiation 
(see equation (14)), M is the mass of the black hole, M is the accretion rate, P = P ra d + -P g is 
the total pressure (P ra d and P g being the radiation pressure and the gas pressure respectively), 
r is the optical depth, v is the kinematical viscosity, k is the Boltzmann constant, m p is 
the mass of the proton, /i is the mean particle mass of the accreting material, a is the 
Stefan-Boltzmann constant, c is the speed of light, / = [1 — (<oGM / c 2 R) l l 2 ] 1, as we are 
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interested in the outer regions of the discs and a is the Shakura-Sunyaev parametrisation of 
the viscosity. Equation (9) gives two different viscosity prescriptions. The first corresponds 
to a situation where a given viscous stress tensor component is proportional to the total 
pressure (discs with this viscosity prescriptions are called a-discs) and the second when 
it is proportional to the gas pressure only (discs with this viscosity prescription are called 
/3-discs) . 

2.1. DISC WITH NO IRRADIATION 

In general we can divide an optically thick disc into four different regimes, namely 

P e < ^rad, K es > K abs inner region (a) 

P g > ^rad, ^es > ^ abs middle region (b) 
P g < -Prad, «es < «abs middle* region (b*) 
Pg > Prad, ^es < ^ abs outer region (c) 

The names of the regions correspond to those used by Shakura & Sunyaev (1973); 
region b* was not included in their paper. We have found that this region exists in discs with 
a^3x 1CT 3 . Which regions are actually present in a given disc depends on M, M, and a. 

Shakura & Sunyaev used only the Rosseland mean for free-free absorption 

Kabs = 6.2 x 10 2 yr~ 3 - 5 cmV 1 . (10) 

However for most relevant densities and temperatures bound-free opacities contribute sig- 
nificantly to the total, and this approximation is quite inadequate. Some consequences of 
using (0) have been pointed out by Cannizzo & Reiff (1992) and Hure et al. (1994b). In 
fact many local considerations are unaffected by the change in opacities, as most quantities 
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depend only very weakly on them (typically as ~ K%bs)- However the global structure of the 
disc is very severely affected by the use of the wrong opacities, as these determine where the 
various regions a - c match to each other. This is illustrated in Figure 1 where the borders 
Rab, Rbc between regions a and b, and between b and c, are given as functions of M for 
two different opacity approximations. The mass of the central black hole is taken equal to 
10 8 M Q and parameter a = 0.001. The solid line for Rb c is obtained using a power law fit to 
the recently compiled solar-abundance opacities (Mazzitelli, 1989 and references therein) in 
the relevant temperature-density regime. 

Kabs = 9xl0 24 pr 3 - 5 cmV- (11) 

The dotted line gives to the R bc values if the opacities are given by equation (10). The filled 
squares are taken from the numerical calculation (Szuszkiewicz, Malkan & Abramowicz 1996) 
where Cox & Stewart (1970) opacity tables were used. The values of R a b are the same in 
both cases, as in regions a and b electron scattering is a dominant source of opacity. We list 
all relevant expressions in Appendix A. 

The disc is terminated at the radius where self-gravity becomes important. This is 
evaluated using the stability criterion for a differentially rotating disc (Toomre 1964) 

Q T = c s Q/nGZ > 1 (12) 

where f2 is the angular velocity. The condition Q T = 1 defines a self-gravity radius 

R sg = (M/vrp) 1 / 3 . (13) 

Knowing the gas density from the given accretion disc model we can find the outer edge of 
the disc. The disc can terminate in any of the regions mentioned above and the appropriate 
density structure must be used. After simple algebra one can find appropriate formula for the 
self-gravity radius for every region in the disc. We give them for all regions and for both a- 
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and /3-discs in Appendix C. We compared our results with the self-gravity radii calculated 
by Hure et al. (1994a) and found reasonable agreement. Many formulae for the self-gravity 
radius used in the literature simply assume that this radius is always located in region a, 
although there is no particular reason to believe this. Region a can also be terminated by the 
change from radiation to gas pressure before the density reaches the self-gravitating value. 
In order to determine in which region, a, b, b* or c a disc actually ends it is sufficient to 
compare the dimension of each region with the appropriate self-gravity radius for given M, 
M and a. The results are shown in Figure 2. 



2.2. IRRADIATED DISCS 

If irradiation dominates the heating, the disc can extend to much larger radii before self- 
gravity starts to be important. A similar case was considered by Ko & Kallman (1991), but 
their results are affected by the use of the opacity form (|To|) . The structure of the irradiated 
disc can be found from equations (0 - |^). Viscous energy generation is now negligible in 
comparison with the flux absorbed from the external radiation, Q\ V v For photons emitted 
radially from a central source, Q nv has the form 

r,Mc 2 (l- X ) D d(H/R) 
Qi " = 4^ R ^R~ (14) 

where 77 (~ 0.1) is the efficiency of the accretion process, x (reasonable values are between 

0.1 to 0.9) is the albedo of the disc and takes into account the projected surface area 

of the disc normal to the radiation flux. If the central emission is not radial but from the 

inner disc Qi rr is multiplied by a factor H/R because this source is foreshortened; cf King, 

Kolb & Szuszkiewicz, (1997). We shall not consider this case here, because we shall find 

that even with the favourable assumption (13) irradiation is never crucial in deciding the 



stability of an AGN disc. We strongly suspect that this is true also for other possible ways 
in which a disc can be illuminated (for example by an extended corona or a jet). However 
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to show this requires additional assumptions about these poorly-understood configurations 
and we leave this for further investigation. For low luminosities and low central masses the 
disc structure is altered by the presence of irradiation. This is shown in Figure 2, where the 
part of the M — M plane for which the disc is terminated in an irradiated region is called 
C + . The thickness of the disc in the irradiated part does not vary with radius as H oc R 9 ^ 8 
(the non-irradiated case) but as H oc i? 45//3S . Irradiation effects are relevant in zone c if 

L M <? _ 7 8 

= ^ 7 x 10- 5 «^M" 5 (1 - X )h (15) 

-^Edd M c 

in zone b if 



^IxlO-VM^fl-^, (16) 



in zone b* (for a-discs) if 



M 3 _5 

^- £ 3 x lO^MKl - X ) J (17) 



and in zone b* (for /3-discs) if 



^-^lxl0V 3 M 8 2 (l-x)o.9 2 (18) 

Here L E dd = 1-5 x 10 46 M 8 erg s -1 is the Eddington luminosity, M c = 2.6 x 10 26 M 8 g s -1 , 
M 8 = M/10 8 M Q and (1 — x)o.9 = 75^- The thin disc approximation we use here requires 
M/M c <^ 0.2. The vertical thickness of the disc in region a is constant and it therefore 
cannot be irradiated by photons emitted radially from a central source, as it is all shadowed 
(cf equation (0)). The situation will be, of course, quite different in the case of irradiation 
by the corona or emission scattered above the disc. 
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3. DISC STABILITY 

Armed with the values of the disc self-gravity radius from the previous sections, we can 
now check the simple criterion T eS (R sg ) > T H for stable disc accretion. We consider values 
Th = 6500, 2000 K, corresponding to the extremes of what is normally claimed for AGN 
discs (Lin & Shields, 1986; Clarke, 1988; Clarke & Shields, 1989; Mineshige & Shields, 1990; 
Cannizzo, 1992). 

Accordingly we identify the stable regimes in regions A, B, B* C, and C + for both a 
and P discs (Figure 2). Here the M — M plane is divided into regions A, B, or B*, C, and 
C + consisting of all models in which the disc is terminated by self-gravity in regions a, b, 
b* , c and c + respectively. The Eddington limit and the limit of validity for the thin disc 
approximation, namely L <^ 0.2Le, are also shown. 

Unlike in the case of low-mass X-ray binaries, even in the extreme case of irradiation 
with very low albedo (90% of the incident radiation absorbed by the disc), irradiation never 
stabilizes the disc. To see this, we note that the criterion for stability in the irradiated disc 
in zone c is 

_ 3 _I 5 M 17 _17 __9_ 24 

7 x 10- 5 a5M 8 2 (1 - X )L Z * 2 - 6c ^ M 8 26 (1 - x)o.f Moo (19) 

IVIq 

where [Th] 20 oo ^ s Th/2000. From this criterion and from the requirement that M/M c <^ 0.2 
it follows easily that irradiation is unable to stabilize the disc in zone c. Equivalent criteria 
for zone b, b* lead to the same conclusion. 

4. DISCUSSION 

Our aim in this study was to investigate whether the thermal-viscous ionization insta- 
bility operates in AGN in the presence of irradiation. We have studied stationary, optically 
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thick, geometrically thin discs, in the range of accretion rates and central black hole masses 
for which these models are self-consistent. It is worth mentioning here that advection domi- 
nated optically thin discs can in principle coexist in some particular regions of the parameter 
space, but which type of the solution will be actually chosen in nature is still an open ques- 
tion. We used a very simple analytic criterion to determine the stability of each model; if 
the disc is hot enough for hydrogen to be completely ionized everywhere all way out till its 
self-gravity radius the ionization instability cannot operate. We identify such hot regions 
in the relevant parts of M - M plane and show them as grey (for a-discs) and hatched 
(for /3-discs) areas in Figure 2. Unlike other authors (Clarke & Shields, 1989; Mineshige & 
Shields, 1990, Cannizzo, 1992) we consider only the upper stable branch of the whole cycle, 
where our method is appropriate. A major advantage of our approach is that we do not need 
a complicated discussion of the limit cycle. This method proved successful in similar studies 
of accretion discs in X-ray binaries. We gave careful consideration to the opacities used 
in our calculations. There are only small differences between results using opacities from 
Mazzitelli (1989) and Cox & Stewart (1970). However differences appear when using simple 
fitting formulae such as (10) instead of (11) (see Figure 1): it is important to check carefully 
that a particular fit found in the literature is appropriate for the range of temperatures and 
densities used in a given problem. 

Another result of our study is that for a ^ 0.003 the region between region a and 
c differs from the standard Shakura-Sunyaev region b. We denote it b*. It is radiation 
pressure dominated, but the main source of opacity is true absorption. We have confirmed 
the existence of this region in numerical calculations of global disc structure performed 
using the Cox & Stewart (1970) opacity tables. It is interesting that b* is stable against 
disc instabilities triggered by radiation pressure (Pringle, 1976): while irradiated it might 
significantly change its properties. 
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In Figure 3 we compare our results with those based on detailed studies of the outburst 
cycle over the parameter space considered by various authors. The dotted lines are from 
Mineshige & Shields (1990), dotted-dashed from Clarke & Shields (1989), long-dashed from 
Cannizzo (1992) and the bold lines from this paper. The short-dashed line gives the Ed- 
dington limit. Our results for non-irradiated disc are in good agreement with those obtained 
previously. Our main result, quite contrary to the case of close binaries, is that irradiation 
does not change the borders between unstable and stable (partially or completely ionized) 
regions. In other words, irradiation by a central point source is unable to stabilize the whole 
disc out to its self-gravity radius. An important reason for this is that one of the effects of 
such irradiation is to move the self-gravity radius even farther out from the central black 
hole. The irradiated disc structure for low-luminosity, low-mass objects differs from that of 
the equivalent discs without irradiation (regions C + in Figure 2). Thus the actual appear- 
ance of the ionization instability might well be affected. This can be studied only by detailed 
calculations of thermal limit cycles in the presence of irradiation. 

We see from Figure 2 that in general AGN discs will be subject to the ionization insta- 
bility, even if they are irradiated by a central point source. For typical AGN luminosities, 
corresponding to central accretion rates ^ lCr 2 M yr _1 , we see from Figure 2 that it is 
inconsistent to assume that the disc is stable. Since central accretion (and thus e.g. X- 
ray emission) is suppressed in the quiescent state, all observed AGN must presumably be 
identified as such only in their outburst states (which last ^ 10 3 yr). Thus AGN currently 
observed to have central accretion rates below the stability limits ~ 1CT 1 — 1O~ 2 M yr -1 
shown in Figure 2 must actually have considerably lower fuelling rates. Even rather brighter 
observed AGN need not be steady systems, but may simply represent the outburst states of 
unstable disc with fuelling rates below the stability limits. 

As pointed out in the Introduction, if most AGN discs are unstable, then in quiescence 
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these systems must be indistinguishable from normal galaxies. Moreover the mass fuelling 
rates needed to power AGN must be much lower than implied by their current luminosities. If 
the duty cycle for the outburst can be made short enough (<; 10~ 2 ), no fuelling rates greater 
than about 10" 2 M yr~ x would be needed in AGN. This would also remove the problem 
that the remnant black holes are predicted to have excessively high masses if accretion is 
continuous (Cavaliere & Padovani 1988). 

Alternatively since most quasars have observed central accretion rates above the sta- 
bility limits in Figure 2, some or all of them could have steady discs. This group would 
then form a separate class with much higher fuelling rates M ~ 0.1 — 1 M Q yr _1 . It is not 
easy to decide between these possibilities by looking at detailed properties of the individ- 
ual systems, as outbursting discs rapidly take on a quasi-steady surface density profile (cf 
Cannizzo, 1993: this property is well known in the context of cataclysmic variables, where 
the persistent systems - novalike variables - look like dwarf novae in permanent outburst). 
A complicating feature is that many of the objects with high steady fuelling rates would be 
subject to the radiation-pressure (Lightman-Eardley) instability. 

We conclude that many (if not all) AGN represent the outburst state of a thermal- 
viscous disc instability. We should then consider candidates for the quiescent state. It is 
tempting to suggest that this may comprise most or all "normal" spirals. Galaxies such 
as our own could therefore harbour moderately massive (10 6 — 1O 8 M ) black holes in their 
nuclei. 
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FIGURE CAPTIONS 



Fig. 1. — The radii R ab) R bc (in terms of the Schwarzschild radius, r g = 2GM/c 2 ) dividing 
regions a and b, and regions b and c, as functions of accretion rate (in terms of the critical 
accretion rate, M c ). R bc shown by the solid line was obtained using opacities given by 
equation (11) and that shown by the dotted line was calculated with opacities given by 
equation (10). R a t are the same in both cases. Here M = 1O 8 M and a = 0.001. The filled 
squares are taken from the numerical calculation (Szuszkiewicz, Malkan & Abramowicz 1996) 
where Cox & Stewart (1970) opacity tables were used. 

Fig. 2. — For a given set of parameters (a, M and M) the global disc structure is uniquely 
determined. Region A contains all models in which self-gravity truncates the disc in a 
region where the pressure is dominated by radiation and the opacity by electron scattering 
(Shakura-Sunyaev region a). Region B contains the models in which the self-gravity radius 
terminates the disc in the zone where the pressure is dominated by gas pressure and opacity 
by electron scattering (Shakura-Sunyaev region b). Region B* is the regime in which the 
self-gravity radius occurs in the zone where radiation pressure dominates, but the main 
source of opacity is absorption rather than scattering. Region C contains models with the 
outer radii determined by self-gravity in the zone where gas pressure dominates and opacity 
is given by absorption (Shakura-Sunyaev region c). In all regions A, B, B* and C the main 
source of heating is due to the viscous energy dissipation. In region C + the main heating 
process is irradiation by the central source instead. The size of this region depends on the 
albedo of the disc, x- It is shown here for two particular values of x- 0-1 an d 0.9. In every 
region the stability properties have been determined, both for a and /9-discs, for two values 
of hydrogen ionization temperature T# (left panel: T# = 6500 K, right panel: T# = 2000 
K) and for different values of a parameter (left panel: a = 10~ 4 , a = 10~ 3 , a = 10~ 2 , 
right panel: a = 10~ 3 , a = 10~ 2 , a = lO^ 1 ). The shaded zones show the location of the 
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a-discs models in which the ionization instability does not operate. The hatched zones on 
the left panel show the location of the /3-disc models in which the ionization instability is 
suppressed. For T# = 2000 K (right panel) the stability properties for a and /3-discs are the 
same, so the shading shows the stable region for both cases. 

Fig. 3. — M - M plane for AGN discs, divided into three distinct regions according to 
the degree of hydrogen ionization: neutral, partially ionized and completely ionized. We 
take a = 0.1. The box with thick solid lines defines the domain discussed in this paper, 
and horizontal line inside it shows the accretion rate above which the discs are completely 
ionized. Previous studies of the M - M plane are shown by different lines: Clarke & Shields 
(1989) - dot-dashed, Mineshige & Shields (1990) - dotted, Cannizzo (1992) - long-dashed. 
The short-dashed line represents the Eddington limit. See text for full details. 
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APPENDIX A 

Size of the regions in non irradiated a and /5-discs 



Here and in Appendices B and C we give expressions for all characteristic radii used in 
our paper. They are given in a form of two equalities. In the first equality, for generality we 
retain a dependence on the electron scattering opacity coefficient k cs = 0.2(1 + X) [cm 2 /g], 
where X is the hydrogen content by mass (in this paper we use X=0.7), and on k q which is 
a constant coefficient used in our fitting procedure: 

«abs = K pT~ 3 - 5 [cm 2 /g] 

The best fit to the solar-abundance opacities (Mazzitelli, 1989) was obtained with k q = 
9 x 10 24 . Moreover, still in the first equality , M should be in [g/s] and M in [g] in order to 
get the radii in [cm]. In the second equality, in order to make the use of our formulae easier, 
we introduced the following quantities: 

M_i = TTT^T-r and M 8 



0.1M C 10 8 M ( 







where M c = 2.6 x 10 26 M 8 [g/s]. All the coefficients in front of these quantities are given in 
[cm]. 

Note that for a less than 0.003 only region b is present in a disc, and for a greater than 
0.003 only region b* is present. 
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R bc = 2.8224 x 10 2 /to M S M§ M§ 



7.9 x KFM^Ms 



i4* c = 1.5089 x 10~ 51 Ko 5 aAMil M§ 



= 4.1 x lO^aTsMi^Mg 15 



18 



R ab = 2.6602 x 10- 17 Kiia^M2T Ma 



1.4 x 10 16 a2iM5iM 8 21 



for a — discs 



_16 10_ n . on -I 

i? a6 * = 3.7322 x 10- 7 /to 45 ftei a^Mis Ms 



= 1.0 x lO^aisM^Mg 45 

for P — discs 



58 _ 20 



R ab * = 4.1353 x 10- 6 «£ko ^a^M^r Ms 



(20) 



o . 36 53 

= 1.0 x 10 16 a-iMi\Mi x 
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APPENDIX B 

Size of the regions c and c + , b and 6 + and 6* and b* + of irradiated discs 



R cc+ = 2.5706 x 10- 28 /to * 5 a^M-Ts Mi 



4 • ± Z2 

= 2.4 x 10 18 a 45 M.^Mg 45 



i? 66+ = 2.0215 x 10- 3 Vs 21 a^M-2T M? 



9 • 

= 4.1 x 10 18 a 5T M.^Mg 21 



for a — discs (21) 



2 ■ 2 , ,13 



_ M 

= 8.8331 x 1Q- 23 K 45 «4%M-f Mi 



2 47 
18^,4- J\>~5 A/f 45 



= 3.6 x lO^a^M.^Mg 4 

for (5 — discs 



R b * b *+ = 3.9392 x 10- 22 Ko ™a&M-% M§ 



= 3.8 x lO^a^M.^Mg 51 
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APPENDIX C 

Self-gravity radius for each region in non-irradiated and irradiated discs 

Region c (for both a and /3-discs) 

2 90 . oo 1 

R sg = 3.9687 x 10 12 K o 15 a45M-45 Ma 

(22) 

22 7 

= 1.8 x 10 17 «lMlfM 8 "« 

Region 6 (for both a and /3-discs) 

R sg = 8.0781 x 10 10 /t! s a^M-^ Ms 

(23) 

= 1.1 x 10 17 «57Mrf M 8 ^ 

Region 6* for a-discs 

^ = 6.5485 x 10~ 12 /tF a§ M^ M5 

(24) 

= 1.0 x 10 17 alM5iM? 

Region 6* for (3— discs 

i? sg = 1.7057 x 10 2 kP a§ M"f M^ 

(25) 

= 1.4 x lO^a^M.^Mg 13 
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Region a for a-discs 

R sg = 4.6733 x KlalMl Ms 

(26) 

= 2.6 x 10 16 a* M^Mg 



Region a for /3-discs 

R sg = 7.6518 x 10- 2 /tlsa3Ms 

(27) 

= 4.5 x lO^asM^Mf 



Region c + 

i? S g = 1.3615 x 10 20 /4aTf M"t m\ 

(28) 

= 1.1 x 10 17 a^MZ['M^ 



Region b + 

R sg = 1.3118 x 10 25 /tia§M^ 

(29) 

= 3.1 x W % aiMZ^MP 
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Region b* + (for a-discs) 

= 8.0208 x 10 14 /4 aiMi M~ 

= 2.1 x 10 13 «- MJ ± M 8 9 

Region 6* + (for /3-discs) 

7.2272 x 10 18 /4 a^M'TE M~ 

1.4 x 10 16 a it MlfM"^ 
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